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Probabilities of a pair of fermions and bosons creation in a static and spatially uniform electric 
field E are represented in the Schwinger formulas by infinite series. It is believed that in weak 
fields the main contribution to the probability is given by the first term of series, however the 
size of the remainder apparently was analyzed by nobody. We study the mathematical structure 
of the Schwinger series by using methods developed during last decades and prove that the given 
series allows an exact summation and the contribution of remainder growths rapidly with the field 
strength. As a rule, it is argued that the pair of particles begin to be produced efficiently from 
the vacuum only in the fields of strength E > E cr . However, the direct calculation shows that the 
Schwinger formula for creation of e + e~ pairs is valid only at the field intensities of E < 0.0291i? cr . 
At higher fields, the probability of pair production in an unit space-time volume exceeds unity. In 
this regard, we refine the formula for the probability of pair creation and numerically find that in 
the field of strength 2.95% of E cr the pair production probability is almost 100%. 



1. More than sixty years ago Julian Schwinger [1] pub- 
lished the formulas that allow to calculate the effects from 
impact of strong electromagnetic fields on the vacuum. 
Were studied the vacuum polarization in strong electric E 
and magnetic H fields, and the resulting one-loop non- 
linear correction to the Maxwell's Lagrange function. Of 
this fundamental paper on quantum electrodynamics we 
are interested in a particular result obtained "inciden- 
tally" by a modest statement of Schwinger |l|. It was 
found that the change in the Lagrangian density causes 
the appearance of imaginary part of the effective action 
that corresponds to a non-zero probability of actual pair 
creation in the field. The probability , per unit time 
and per unit volume, that a pair of particles carrying spin 
s is created by_a homogeneous and static electric field is 
expressed by 
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where 1 a = e 2 /An, E cr — m 2 /e is the so-called critical 
field, m the particle mass, e the charge on the particle. 
The '+' sign corresponds to creation of fermion pairs, 
and the '-' sign for production of bosons. At production 
of e + e~ pairs, E cr = 1.32 x 10 16 V/cm. The scale factor 
E cr was introduced and evaluated for the first time in 
1931 by Fritz Sauter Q. Equation (fTJ) corresponds to 
the weak field approximation. 

The peak electric field achieved at the focus of today's 
most powerful lasers still several orders of magnitude 
lower E cr . This is the main obstacle for a direct ex- 
perimental verification of given fundamental mechanism 
of particle production. Despite of this, the Schwinger 
mechanism from QED, was the basis for successful phc- 
nomenological models of hadron production in collisions 



at high energies via the tunnel creation of gg-pairs in 
string-like chromoelectric field 

From the form in which presented the result ([T]) and 
the value of E cr , it seems self-evident that in weak fields 
the main contribution to the production probability gives 
the first term of the series 
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and this probability is very small. It should be noted, 
however, that the total value of residual terms of an infi- 
nite series, apparently not analyzed. One of the purposes 
of this article 2 is an evaluation how strongly change the 
remainder of the series (JIJ with increasing field strength 
. The analysis of this formal mathematical problem al- 
lowed to find the exact sum of the infinite series ([TJ, 
and find a number of unexpected consequences, not men- 
tioned by other authors. 

2. To solve the above problem we introduce the re- 
duced field /3 = E/E cr , and denote the exponential fac- 
tor in (pQ) by x = exp(— n/f3). As a result, (UJ takes the 
form 
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Analysis of the mathematical literature has shown that 
the series in ([3]) belongs to a class of polylogarithms, and 
in our particular case we are dealing with the diloga- 
rithm. In mathematics, the dilogarithm studied since 
Leibniz and Euler, but in the last quarter of the 20th cen- 
tury polylogarithms and hyperlogarithm again attracted 
attention and appeared in many branches of mathematics 
and physics. Polylogarithms and reference information of 
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TABLE I: The exact values of the Rogers dilogarithm 
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them are included in the standard set of special functions 
in most packages of applied mathematical programs. 

Let us write out the definition of the dilogarithm and 
its integral representation following the interesting orig- 
inal review [fj, 
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Thus, in terms of the function Li 2 (x) the probability of 
pair creation we rewrite as follows, 
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and the problem of summing the infinite series in ([T]) can 
be regarded as solved. However, to clarify the meaning of 
this result, let us recall some properties of the dilogarithm 
[f|. In this case it is more convenient and more compact 
to use the Rogers dilogarithm L(x), < x < 1, defined 
as 



L(x) = U 2 (x) + ~ ln(x) ln(l 



x). 
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The function L(x) allows an analytical continuation into 
the complex plane with cuts (— oo,0] and [1, +oo) along 
the real axis. Using a number of theorems and functional 
relations, L(x) can be extended to the whole real axis. 
There are only few real points {xi}, where the values of 
L(x) are known exactly. Their values are summarized 
in Table I. In general, for more than two centuries of 
research were opened a large number of functional iden- 
tities, allowing to relate dilogarithms with different pow- 
ers of x. There are identities and for more complicated 
rational arguments. 

Coming back to the problem of calculating the prob- 
ability of pair creation by (0, we summarize that by 
using the integral representation of Li2(x), one can cal- 
culate tJ/ s )(/3) with any given accuracy, controlled by the 
precise values of L(x^) from Table I. Another approach is 
possible also. Having determined numerically L(x) for a 
given Xo, the values of L(x) at other points are restored 
through a chain of functional identities. 

To conclude this section we estimate, at which values 
of j3 the remainder of the series (JXJ) starts to give a signifi- 
cant increment to ©. We consider only the case of e + e~ 
pairs. To do this, compose a ratio of the total Schwinger 




FIG. 1: e + e~ pair production. The ratio of the total 
Schwinger series {!]) to its first term |2| as a function of the 
reduced field j3. 



series to its first term 
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By a numerical integration one find (see Fig. 1) that only 
at /3 > 0.5 the remainder of the series begins to play an 
important role, making the correction of 6% at f3 = 2. 
But these values of /3 are beyond the weak field approx- 
imation. Thus, ^) indeed is a very good approximation 
to (J). 

3. The method of the simplified record of results when 
a group of some physical constants are set equal to unity 
is widespread in theoretical physics. In the case under 
study ([J) it is % = c = 1 . As noted by several authors [3] , 
"simplifications" achieved in the natural system of units 
(the system of Planck and Hartree, a relativistic system), 
sometimes are fraught with illusions 3 and Schwinger for- 



3 In this regard, in the book Q on page 97 there is a remarkable 
comment: "Thus, simplification of the formulation of physical 
laws by using the natural systems of units is bought at the cost 
of replacing the physical equations by numerical one. In the used 
units, the numerical equations relates among themselves not the 
physical quantities but only their numerical values. A possibil- 
ity of comparison and check of the dimension of the considered 
physical quantities is lost. Due to preservation of the same let- 
tering a substitution of physical equations by the numerical one 
is masked and there is an illusion of operating with real physical 
quantities and equations. Arises abstractness of used equations 
and the theory as a whole." 



3 



mula (|2|) may serve as a vivid illustration. 

In the vast literature devoted to pair production in 
electromagnetic fields, we have not ever met a graphical 
representation of ([2]). Before we draw the graph is' 1 / 2 ) 
as a function of /3, one need to choose the system of 
units. In [l[ Schwinger likely used the Heaviside system 
of units, the most common in scientific literature of the 
time. In @ the dimension of the probability density 
is cm _3 -s _1 and the dimension of electric field V/cm. 
We need to restore the "invisible" conversion factor to 
match the dimensions on both sides. Let us reproduce 
schematically the logic of calculations led to ([1]) in order 
to recover all constants h and c and bring the final result 
in the SI units 4 . 

Suppose we observe during a period of time AT for 
a three dimensional space volume AV^ filled with a ho- 
mogeneous and static electric field E. Wc write the 
wave function of the vacuum as $ = exp(iS/h), where 
the action S is related to the density of the Lagrangian 
function L via S — C AVAT. As a result, the proba- 
bility that the vacuum will be stable is W v = l^l 2 = 
cxp(— 2lm£ AVAT/h), and the probability of pair cre- 
ation in the field is 



W p = 1- W v = l-exp(--Im£AVAT). (8) 

n 

In the weak-field approximation, keeping only the first 
term in the expansion of the exponential, we reproduce 
the Schwinger result, 
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that is, the probability of pair production per unit time 
per unit volume is w^ s ' — 2lm C/K. The dimensions of 
the functions S and C are fixed by the above formulas 
allowing to restore the required number of the desired 
constants. Thus, E cr = m 2 c 3 /eh and 
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where 
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For e + e~ pairs, after substituting numerical values of 
the constants one find that in ([2]) behind the seemingly 
small pre-exponential factor is "hidden" the very large 
dimensional constant wq = 1.087036 x 10 50 cm -3, s -1 . 
Consequently, the earlier qualitative conclusion about the 
smallness of w' 1 / 2 ' may not be entirely correct 5 . 



In Fig. 2 by open dots shown Wp 1 ^ 2 ^ calculated with 
© and (JTOj) with ATAV = 1 cm 3 -s. Taking into account 
the magnitude of Wq, the result is not so surprising, de- 
spite the fact that the field strength is much lower than 
the critical value. As follows from Fig. 2, the Schwinger 
formula is valid only for fields of E < (3 c E cr . At higher 
field, the probability of pair production exceeds unity. 



From (fT0|) we find /3 C at which Wp 
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Let us denote A = ? \/wqATAV and sequentially making 
change of variables y/woATAV (3 C = z and z 
get an equation 



W(A)e w ^ = A 
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for the Lambert W(A) function. The Lambert function 
has the same fate as that of the dilogarithm. Interest in 
it was revived in the last quarter of the 20th century [l(| , 
ll| . The required (3 C is obtained from 
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For numerical estimates we use the asymptotic expansion 
of W(A) 
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where L x = In A, L 2 = In In A. Thus, p c (W p = 1) « 
0.02905 at AVAT = 1 cm 3 -s. 

The reason on which at j3 > f3 c the unitarity is violated, 
lies in the use of the weak field approximation leading to 
©. To restore the unitarity, it is necessary to take a step 
back to ©, which is free of the problem. In ([5]), 2Im£/?i 
should be replaced by from ([5]) and take into account 
the new system of units (ITU|) . After completing the nec- 
essary substitutions, we obtain the corrected formula for 
the probability of pair creation in the volume AV 



WW(J3) = 1 - exp{±(s + -)u, /? 2 Li 2 (±e-^) A^AT} 

during the observation period AT. In (|14[) already makes 
sense to substitute the exact form we found for the 
Schwinger series, because we are not limited by the weak- 
ness of the field 6 . It should be emphasized that (|14|) ex- 
plicitly demonstrates how the probability of the process 
depends on the size of the space-time volume. But this 



4 However, in the following, as the unit of length is more convenient 
to use 1 cm. 

5 Our expression for wq coincides with the analogous scaling con- 
stant of Refs mm. 



6 Provided that the two-loop correction [12| | to the Heisenberg- 
Euler Lagrangian density gives a small contribution. 
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FIG. 2: The probability of an e + e~ pair production in an unit 
4-volume, AVAT = 1 cm 3 -s, as a function of the reduced 
field p. The open dots shows the result of calculation by 
equation @, violating the unitarity. The solid curve shows 
the calculation by (fl4j) . 



dependence is extremely weak if compared with the de- 
pendence on E. For example, if at a given field strength 
the probability of pair creation is 95%, and we wish to 
reduce the vacuum field by the factor ten, keeping the 
same particle production rate, then it is necessary to in- 
crease the field space-time volume in 10 2 -exp(^) w 10 420 
times, that is, literally to the cosmological size. 

In Fig. 2 by the solid curve shown the probability of an 



e + e pair production calculated by means of (IT4|) with 

AVAT — 1 cm 3 -s. As is evident from the figure, Wp 1 ^ 2 ^ 
increases rapidly within a very narrow range of variation 
of E. Indeed, when the field strength is 2.75% of E cr the 
pair production probability is less than 1%, but in the 
field of strength 2.95% of E cr , e + e~ pairs are produced 
with almost 100% probability. This behavior is typical 
for tunneling processes. 

4. In the present article we prove that the Schwinger 
series ([I]) for the probability of a pair of bosons and 
fcrmions production in a uniform and static electric field 
allows to be represented in terms of the dilogarithm. In 
a weak field, j3 < 0.5, this series with a very high ac- 
curacy is described by the first term. With increasing 
/3 the remainder of the series begins to play an impor- 
tant role, making the correction of 6% at j3 = 2. The 
restoration of all constants Ti and c in the Schwinger for- 
mula reveals a very large pre-exponential factor wq = 
1.087036 x 10 50 cm~ 3 -s _1 . Its presence determines a 
sharp increase of the pair production probability within 
a very narrow range of field intesities and in the weak 
field approximation leads to a violation of unitarity at 
E > 0.0291EW. Equation (fi"4"f corrects this defect, sat- 
isfies the unitarity condition and includes the space-time 
volume of the process. Numerical estimates show that 
in the field of strength 2.95% of E cr within the four- 
dimensional volume of 1cm 3 -sec, the probability of an 
e + e~ pair production is close to 100%. 
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